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Abstract 

We examine contributions from Majorana phases to lepton flavor violating processes in the 
framework of the minimal supersymmetric standard model with heavy right-handed neutrinos. 
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t^J" , All phases in the complex neutrino Yukawa matrix are taken into account in our study. We 

CN ' find that in the scenario with universal soft-breaking terms sizable phase effects can appear on 

branching ratio of u — > can be considerably enhanced due to the Majorana phases, so that 
Q_ij it can be much greater than that of r — > ^7. 
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1. INTRODUCTION 



In the standard model lepton flavor violation (LFV) is negligible, while it can be 
sizable in new physics models such as those based on supersymmetry (SUSY). Therefore 
search for LFV can be a good probe of new physics. Observed tiny neutrino masses 
may be explained by the seesaw mechanism^], assuming heavy right-handed Majorana 
neutrinos, which are compatible with the scenario of grand unified theories (GUTs). In 
the framework of SUSY models, LFV is induced through one-loop diagrams with slepton 
mixing (2)]. In the SUSY model with right-handed neutrinos, the slepton mixing can 
be induced from the renormalization group effect of the neutrino Yukawa interaction 
between the scale of right-handed neutrino masses and the GUT scale, even when soft- 
SUSY-breaking terms are universal at the GUT scale. 

The neutrino mass matrix obtained via the seesaw mechanism generally includes two 
Majorana phasesj^j. They can be directly searched through neutrinoless double beta 
decays jj]. The existence of these Majorana phases can play an important role in various 
phenomena such as leptogenesisjsj], lepton number violating processes and so on. Searches 
for these phenomena could provide a hint for the neutrino Majorana mass matrix. Fur- 
thermore, as we shall show below, the prediction on LFV can be drastically changed by 
the Majorana phases. 

In the present paper, we explore LFV processes such as fi — > in the framework 
of the minimal supersymmetric standard model with right-handed Majorana neutrinos 
(MSSMRN) under the assumption of universal soft-SUSY-breaking terms at the GUT 
scale Mqut- Neutrino mass matrix m v is given by m v = Yj D^Y^l) 2 , where Y v is the 
neutrino Yukawa matrix, Dr is the right-handed neutrino mass matrix which is diagonal, 
and <p° u is the neutral component of the Higgs doublet with hypercharge —1/2. In the 
basis where the charged lepton mass matrix is diagonal, the neutrino Dirac mass matrix 
mo = Yv{4>V) can be parameterized by[f], Q] 



DbRx/LLU^ 



where _D„ is the eigenmatrix of neutrino masses, R is a complex orthogonal matrix 



(R R = RR = 1), and U is the neutrino mixing matrix. In Refs. 



the decay 



rates of li — *> (i ^ j) are evaluated by assuming that R is a real orthogonal matrix 
and that the right-handed neutrino masses are degenerate; i.e., Dr — M x 1 where M 
is the heavy Majorana mass scale. Under this assumption, the effect of Majorana phases 
on the low energy phenomena is screened. The relation among the branching ratios is 
given by 

Br(/x - erf) ~ Br ( r Z? ' ^ « Br (r -> py) , (2) 
Br(r — > v e v T e) 

where current neutrino data have been used. The hierarchical Dn case with a real R has 
been a— in Re, Q. On the other hand, the lm p„ rt an C e of the t—t of « as a 
complex matrix has been pointed out in Ref. [3], by showing that phases in R can give a 
substantial effect on low energy phenomena. 

In this paper, we discuss the role of the imaginary part of R, and study the combined 
effect with Majorana phases in neutrino mixing matrix on the branching ratios of the 
LFV processes. We assume that Dr = M x 1. We obtain analytic expressions of the 
branching ratios in two limiting cases: i.e., one is the case with R being approximately a 
real orthogonal matrix, and the other is with R being a typical complex matrix. We find 
that 

Br(/x - erf) ~ ^ - » Br(r - /i 7 ) , (3) 

in the wide range of the parameter space for a typical complex matrix R. The branching 
ratio of fi — > e 7 can be enhanced in comparison with that of r — >• ^7. This is a novel 
feature with a complex i?. We also give numerical calculations in order to see how these 
two limiting cases are extrapolated. 



2. EVALUATION OF LFV BRANCHING RATIOS 

In this section, we briefly review LFV in the MSSMRN, and discuss the Majorana 
phase effects on LFV processes. 

In the model based on SUSY, LFV processes can occur at the low energy scale through 
the slepton mixing. In the MSSMRN, sizable off-diagonal elements of the slepton mass 
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matrix can be induced by renormalization group effects due to the neutrino Yukawa 
interaction between Mqut and M, even when universal soft-breaking masses are assumed 
at Mqut- The induced off-diagonal elements are approximately expressed as[" 



K) 



6m 2 , + \ A o I 



In 



Mi 



GUT 



(4) 



16vr 2 ~ M 

where mo and A are universal soft-SUSY-breaking parameters. The decay rates for LFV 
processes £{ — > (i ^ j) are given by 



m - in) 



a 



EM m t 



IKkl 2 



192tt 3 m 



tan 2 j3 



(5) 



where «em is the fine structure constant, rn SUSY represents the typical mass scale of 
SUSY particles, and tan/5 is the ratio of vacuum expectation values of the two Higgs 
doublets. The branching ratios are related to each other as 



Br( / u — > ej) 



i 



{m ] D m D ) l2 \ 2 



Br(r 
Br(r 



Br t 



\{rn} D m D ) 2 z\ 



5.6 x 



\{rrJ D m D ) 



12 



{m} D m D ) 2 3\ 2 



ej) \{m) D m D ) lz \ 2 



(6) 



0.1784 is used. These ratios are determined 



Br(r ->// 7 ) |(m^m D ) 23 | 2 ' 

where experimental result Br(r — > v e v T e / 
only by the neutrino Yukawa matrix. 

We work on the basis that the right-handed neutrino mass matrix is diagonal, and 
assume that the matrix is approximately proportional to the identity matrix; i.e. Dr ~ 
M X 1. By using Eq. JTJ) we obtain 



m} D m D ~ MU ^/D V R ] R^/D y U ] = MU ^D V Q ] 'Qy/D^lP . 



(7) 



Here we have introduced a real orthogonal matrix O by R = OQ, where Q is a product 
of Q a (a = 1-3) with 



cosh yx i sinh yi (A 



Qi 



-i sinh yi cosh y\ 
1 



(l 



Mi 



o 







cosh y 2 i sinh y 2 
^0 — zsinh?/2 cosh?/ 2 y 



g 3 



cosh ys i sinh 2/3^ 
1 
^— isinh?/3 cosh?/3y 
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The matrices Q a satisfy that Q\ = Q a and Q 2 a (y a ) = Qai^Ua)- The matrix Q plays a 
role not only to introduce the complex phases but also to change the size of Yukawa 
couplings.* 

The neutrino mixing matrix U is separated into two parts, U = Umn &IL where Umns 
is the Maki-Nakagawa-Sakata matrixj9( in the phase convention of Ref. [10j and P is the 
Majorana phase matrix given by P = diag(l, e ia °, e 1 ^ ) with «o and (3q being Majorana 
CP violation phases JjJ . In order to see qualitative features, we here take the Bi-maximal 
mixing solution ll| 



/ i 



j jBi—max 
MNS 



V2 v/2 
_ 1 1 1 

2 2 ^2 
1 1 1 

V 2 2 V2/ 



(9) 



for analytic calculations. In particular, we consider the following three cases for D u ; 
the normal hierarchical (NH) case (mi <C 7772 <C m.3), the inverse hierarchical (IH) case 
(rrt3 <mi~ m, 2 ), and the quasi-degenerate (QD) case (mi ~ 777-2 ~ ^3); 



NH : mi ~ , m 2 — \ Am| , 777 3 ~ a/ Am 



atm 



QD : 777i = 777 , ?77 2 ~ 771 + 



Am| 



m 3 ~ m + 



Am 



atm 



Here, Am| = m\ — m\(= 8.0 x 10 5 eV 



(10) 



(1 Aw 2 
1 ~ 2 Am 2 ° ' ' m2 - V 7 Am2 tm , m 3 ~ , 



;i2) 



2m 2m 

■ is the squared mass difference for the solar 



neutrino mixing, and Amf tm = \m^ — m%\(= 2.5 x 10 eV is that for the atmospheric 
neutrino mixing. 

To evaluate rn : D m D , we consider the following two limiting cases. 



(a) The small y a limit (R is real.) : 
We have Q = 1, and thus m^m^, 



Mu mns ^DvUmn™ aX \ where the elements of 



We note that fine tuning of order 0(e Va ) is necessary to obtain the light neutrino mass scale in the 
case of y a > 1. 



5 



m D m D are determined by neutrino masses and the mixing matrix as 



{m ] D m D ) 12 = {m ] D m D ) lz 
M 



M 



2V2 



(m 2 - mi) 



{m^ D m D ) 2 3 = —(mi + m 2 - 2m 3 ) . 



(13) 



We then obtain from Eq. JBJ) that 

f 5.6 xff^- 



Br(/x — > e7) 
Br(r -> /i 7 ) 



5.6 x § 
5.6x±( 



Ami N 2 



Am atm 

Ami 



~0.23 for NH 
7.7 x 10~ 4 for IH 
3.1 x 1(T 3 for QD 



(14) 



For all the cases, it turns out that Br(/x — > e'y) ~ 5.6 x Br(r — > ej) <C Br(r — > ^7), 
as pointed out in Refs. [3, B| - in this limit, the Majorana phases do not affect the 
LFV processes. 

(b) The large y a case : 

The matrix Q has a simple form. First, the matrices Q a behave as 



Ala 



(15) 



where 



(1 ,o^ 



Qi 



V2 



-i 1 




A) (A 



Q 2 



V2 



1 % 
-i 1 



^ 1 z * 



v/2 




-z 1 



(16) 



They satisfy Q\ = Q a and Q 2 a = s/2Q a . As for the product of Q a such as Q E 
{QaiQbQa-, QcQbQa}, we find an interesting relation as 

Q j Q = V2Q a . 



(17) 



By using Eq. fTT]) QtQ is expressed by 

g2(j/l+J/2+?/3) 



4V2 



Qa- 



This means that Q^Q is characterized by three independent matrices Q a {a = 1- 
3) for large ?/ a . Thus, we examine the following three cases, taking R = OQ a ~ 
e*>OQa/V2. 



(b-1) R = OQ 1 
We have 



{m ] D m D ) l2 = -{m j D m D ) 13 



f Me 2yi \ m 2 — mi + i1^Jm\m 2 cos a>Q 
\~2~J 271 



{rrJ D m D ) 23 



/ Me 2yi \ mi + m 2 — 2 v /mim 2 sin a 

I ~2~~ J 4 • (19) 



Thus the LFV branching ratios are related as Br(/x — > e7) ~ 5.6 x Br(r — > ej) 
for all cases. For the NH case, we obtain Br(/x — > &y) — 11.2 x Br(r — > /17). 
For the IH and the QD cases, one finds 

Br ^ - 67) x - 11.2 x . (20 ) 

Br(r — > pij) (l + sina ) 

This ratio is a function of ocq. It is 11.2 for a = or tt, and for «o = 7r/2. 
(b-2) i? = OQ 2 

The difference of the Majorana phases «o — ft enters into 7n i D 7n D , 

. t . f Me 2y2 \m 2 + iV2W7n 2 m 3 e^ ao - f3o) 
(m D m D ) 12 = _ , 

\^~) 271 ' 

(Meto\ -m 2 + 2m 3 + 2v / 2z v ^ 2 T^cos(a -/?o) , 01 . 
{m D m D ) 23 - I — - — I . [21) 

For the NH case and the IH case, the branching ratios of i{ — *> £j~f are related 
to each other as Br(/i — > e'-f) ~ 5.6 x Br(r — > e7), and 



(m^m D ) 13 



Br(/z - e 7 ) I 5.6 x \l 4^r- ^ 1.0 for NH 



(22) 



Br(r ^/i 7 ) " n 2 for ffl 

For the QD case, we obtain 

Br(/x -> e 7 ) 3 + 2 v/2 sin(a - ft) 

Br(r — > ^7) 1 + 8 cos 2 (ao — ft) 

Br(r -> e 7 ) _ 2(3 - 2^ sin(a - ft)) 

Br(r — > /ry) 1 + 8 cos 2 (ao — ft) 
We have Br(/i — > cy) ~ 5.6 x Br(r — > e7) ~ 3.7 x Br(r — > 7x7) for a — ft = 
or 7T. The ratio Br(/i — > e7)/Br(r — > /ry) takes its minimum value 1.9 at 
"o - ft ^ -vr/2. 

7 



(23) 



Br(/x->e7) 
Br(r^/ry) 



Br(r— *erf) 



NH 



IH 



QD 



NH 



IH 



QD 



small y a 
Q = l 



0.23 



7.7 x 1(T 4 



3.1 x 10" 



5.6 



Qi 



11.2 



11.2 x " os a \. 



5.6 



large y a 

Q 2 



1.0 



11.2 



1 1 o v 3-2v / 2sin(ao-/3o) 
1J - Z * l+8cos^(a -/3 ) 



5.6 



5.6 x 



3-2-^/2 sin(ao-^o) 
3+2 N /2sin(a u -^u) 



Q 3 
< l 



11.2 



ii 2 v 3+2\/2smJo 
1+8 cos^ /3q 



5.6 



°' U 3-2^8111/30 



TABLE I: Summary of the ratios of the LFV processes. 

(b-3) i2 = OQ 3 

The Majorana phase (3q enters into m} D m D . We obtain 

f Me 2y3 \ -mi + i\pl^m x m z e- i(io 



{m' D m D ) l2 
{m ] D m D ) 



2V2 



■iPo 



13 



m D m Z))23 



v 2 

/Me 2y3 \ mi 
1^7 _ 2v/2 
f Me 2m \ -mi + 2m 3 - i2V2 y /m 1 m 3 cos /3 



(24) 



In this case 



mt,?i 



Diib D )ij can be obtained from case (b-2) by replacing m2 
with mi and «o ~~ A) with ir — (3q. The branching ratio Br(yU — > e7) is 
suppressed in the NH case because of Br(/i — > e^y) ~ 5.6 x Br(r — > e^y) ~ 
5.6 x (mj/Am^ tm )Br(r — > ^7) <C Br(r — > ^7). For the IH case, the branch- 
ing ratios are related to each other as Br(/i — > ej) ~ 5.6 x Br(r — > e^) ~ 
11.2 x Br(r — > ^7). For the QD case, relation among the ratios of branching 
ratios is obtained from Eq. (j2*3j) by changing a — Po to n — /3 . 

The results are summarized in Table 1. 

For the ratio Br(/i — > e7)/Br(r — > ej), Q does not contribute except for the QD case 
with Q = Q 2 or Q 3 where the Majorana phases give a significant effect. The drastic effect 
occurs for Br(/i — > &-/) /Br(r — > /17) by Q or by the interplay between Q and the Majorana 



S 



phases. The substantial enhancement arises in Br(/x — > ej), which is a quite different 
feature from the case with Q = 1. By the introduction of Q(^ 1), the Majorana phases 
can affect the LFV processes. This fact is a quite interesting because the observation of 
the LFV processes would give useful information of Majorana phases. 

3. NUMERICAL RESULTS 

In the previous section, we consider the two limiting cases for the parameter y a . When 
y a take the intermediate values, we may guess the result by extrapolating from the two 
limits, but some non-trivial structure might appear. Therefore, we perform the numerical 
evaluation of the LFV branching ratios for three typical cases, R = OQ a (a = 1-3). 
Neutrino mixing parameters are taken to be tan 2 6* = 0.45^^, sin2# aim = and 
sin# 13 = 0. The values for M and Mqut are taken as M = 10 10 GeV and Mqut — 
2 x 10 16 GeV. The SUSY parameters are taken to be m = A = m SUSY = lOOGeV 
and tan/3 = 10. For standard model parameters a EM = 1/137 and v = 246GeV are 
used. It will be shown that the ratios of the branching ratios are not sensitive to SUSY 
parameters, right-handed neutrino mass scale, and the GUT scale. 

We analyze the y a dependences of Br(/x — * ej) /Br(r — > fi-f). The result for the NH case 
is shown in Fig. 1. We find the smooth extrapolation in Br(yU — > ey) /Br(r — > ^7) between 
0(0.1) and 0(1) for R = OQ 2 with a - A> = and between 0(0.1) to O(10~ 6 ) < 1 
for R = OQ 3 . For R = OQi, some structure is observed between O(0. 1) and 0(10). 
The ratio blows up around y ~ 1.3 due to the vanishing Br(r — > fiy). There is no a 
dependence. 

The IH case is shown in Fig. 2. The dotted (dashed) curve represents Br(/i — > 
e7)/Br(r — > fiy) for R = OQ2 (OQ3) where the smooth extrapolation is found between 
very small value to about 50 (2), where there is no «o dependence. The case R = OQi is 
shown for solid curves, which has the Majorana phase ao dependence. For all cases, we 
find the smooth extrapolations between two limiting values, the small y a and the large 

Va- 

For the QD case with R = OQ2, the ratio of the branching ratios depends on «o — A>, 
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FIG. 1: The ratio of the branching ratios is shown in the NH case for Q± (solid curve), for Q2 
with oq — 0o = (dotted curve) , and for Q3 (dashed curve) . 




0.5 1 1.5 2 2.5 3 



Va 

FIG. 2: The ratio of the branching ratios is shown in the IH case for Q\ with a$ = 0, 37r/4, 2>n/2 
(solid curve), for Q2 (dotted curve), and for Q3 (dashed curve). 

and is roughly obtained by replacing (3q to 7r — (olq — (3q) in the formula for R = OQ%. 
The results for R = OQ\ are similar to those for the IH case with R = OQ\. In Fig. 3, 
we show the y 3 dependence for the case with R = OQ 3 for a — (3 = 0, 37r/4, 37r/2. The 
enhancement occurs for a — (5q = 37r/4 because Br(r — > fij) is suppressed. 

In Fig. 4, Br(/i — > e^y) with R = OQ\ is shown as a function of y±. As yi grows, the 
neutrino Yukawa couplings become large for all the neutrino mass spectrum. Thus, the 
smooth extrapolation is obtained, so that the two limiting cases give the general trend 
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00 = 3tt/2 ; 

A) =_o_ _ " 

A) = 3~7r/4~ ■ 

0.5 1 1.5 2 2.5 3 

2/3 

FIG. 3: The ratio of the branching ratios is shown in the QD case for Q% with (3q = 0, 37r/4, 2>tt/2. 

l.xlO" 9 
1.x 10" 11 

T 

a. l.xio -13 

! — i 

1. x 10" 15 

0.5 1 1.5 2 2.5 3 

FIG. 4: The LFV branching ratios Br(/i — ► e7) are shown in the NH case, the IH case with 
«o = 37r/4 and the QD case with «o = 37r/4 for R = OQ\. 

of the y a dependence. In many cases, the Majorana phases play an important role on 
the prediction of the LFV processes. Therefore, we can obtain useful information of the 
Majorana phases from the experimental data of the LFV processes. 

4. CONCLUSION 

We have shown the importance of the complex nature of the neutrino Yukawa matrix 
for the case of the degenerate right-handed neutrino masses. With the complex R, the 
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Majorana phases play an important role for the prediction of the LFV processes. In 
order to see the effect analytically, we have taken the parameterization, R = OQ. We 
have considered the two limiting cases; the small y a case with Q = 1 and the large y a 
case with complex matrix Q. We have obtained the analytic expressions for ratios of the 
branching ratios of fi — > ej, r — > /ry and r — > ej, which are shown in Table 1. The effect 
of Q is sizable and gives enhancement of Br(/i — > e7)/Br(r — > ^7) in many cases. In 
particular, the Majorana phases contribute to some cases. This would give a possibility 
to obtain useful information of Majorana phases by observing the LFV processes. This 
is quite interesting and important because extracting the information for the Majorana 
phases can be used to examine the nature of neutrinos. 

It may also be interesting to discuss the possibility to determine the neutrino Yukawa 
matrix by analysing the double beta decay, the \i~ —>■ e + jj, ^| and fi~ — > fi + conversion 
|l5l | . the LFV processes which occur through SUSY contributions, and the leptogenesis. 
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